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O ■ ABSTRACT 
(N 

\ Context. Although deuterium enrichment of water may provide an essential piece of information in the understanding of the formation 

. of comets and protoplanetary systems, only a few studies up to now have aimed at deriving the HDO/H2O ratio in low-mass star 

forming regions. Previous studies of the molecular deuteration toward the solar-type class protostar, IRAS 16293-2422, have shown 
that the D/H ratio of water is significantly lower than other grain-surface-formed molecules. It is not clear if this property is general 
or particular to this source. 

Aims. In order to see if the results toward IRAS 16293-2422 are particular, we aimed at studying water deuterium fractionation in a 
second low-mass solar-type protostar, NGC1333-IRAS2A. 
' Methods. Using the 1-D radiative transfer code RATRAN, we analyzed five HDO transitions observed with the IRAM 30m, JCMT, 

I and APEX telescopes. We assumed that the abundance profile of HDO in the envelope is a step function, with two different values in 

\^ , the inner warm (T>100 K) and outer cold (T<100 K) regions of the protostellar envelope. 

Results. The inner and outer abundance of HDO is found to be well constrained at the 3cr level. The obtained HDO inner and outer 
fractional abundances are x™° = 6.6 x 10"**-1.0 x 10"'(3cr) and x™° = 9 x 10""-1.8 x 10"'(3cr). These values are close to those in 
IRAS 16293-2422, which suggests that HDO may be formed by the same mechanisms in these two solar-type protostars. Taking into 
^ , account the (rather poorly constrained) H2O abundance profile deduced from Herschel observations, the derived HDO/H2O in the 

$_( ' inner envelope is > 1 % and in the outer envelope it is 0.9%-l 8%. These values are more than one order of magnitude higher than what 

^ ' is measured in comets. If the same ratios apply to the protosolar nebula, this would imply that there is some efficient reprocessing of 

I the material between the protostellar and cometary phases. 

Conclusions. The HtO inner fractional abundance could be further constrained by an analysis of newer observations of high-energy 
H'^O lines. This new observations would be required to understand water fractionation in more detail. 

^ ■ Key words. Astrochemistry - Stars: solar-type - Stars: formation - ISM: molecules - ISM: individual objects: NGC1333-IRAS2A 
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2^ I 1. Introduction 2422. This can be interpreted in terms of the chemical for- 

" ■ . mation pathways of these molecules - both are believed to 

^ Molecular deuteration has been the focus of many studies, ^e grain surface products in cold dark clouds thro ugh the for- 

O The abundance of deuterated molecules has been shown to be nation process of successive CO hydrogenation (iParise et al.l 

^ : significantly increased in low-temperature environments, com- Because water is the main constituent of grain mandes, 

^ . pared to the D/H elemental ratio in the interstellar medium st^^ies have been performed to investigate if deuterium enrich- 

■ (~1.5xlO-MLinsy 120031). Over the last decade, very high j^ent of water is similar to the fractionation of formaldehyde 

^^^^""^^ molecules have been found in low-mass ^nd methanol. Several HDO transitions have been observed to- 

p><; protostars, e.g. D.CO/H^CO^O.OS jCeccarelh et al.||1998|) and ^^^^ jr^s 16293-2422, and were modeled to derive th e HDO 

^ ; CD3OH/CH3OH=0.01 ( |Parise et al. || 2004|) toward IRAS 16293- abundance profile in its envelope (iParise et al.l 12005^ . Using 

■ ■ ■ ' the H2O abun dance profile as derived from ISO observations 

* Based on observations with the APEX telescope, the IRAM (ICeccarelh et al. 2000). the HDO/H2O ratio was found to be 

30m telescope and with the James Clerk Maxwell Telescope ~ 0.03 in the inner hot corin o, and <0.002 in the cold outer 

(JCMT). APEX is a collaboration between the Max-Planck-Institut envelope (IParise et al.ll20o!ah . i.e. significandy lower than the 

fiir Radioastronomie, the European Southern Observatory, and the formaldehyde and methanol fractionation in the same source. A 

Onsala Space Observatory. IRAM IS a European collaboration between , ^ tit-,/-,- • 1 -.u 1 f 

. ^xTTir. XT ■ , , , , , ■ -n X , ,Tr^^ search for HDO ices in low-mass sources with large Iractiona- 

the CNRS (Centre National de la Recherche Scientifique), the MPG . „„ ,,,, , , , ,,t^^ • 

(Max-Planck-Gesellschaft) and the Spanish IGN (Instituto Geogrfico ^lon of formaldehyde and methanol found no HDO ices to a very 

Nacional). The JCMT is operated by The Joint Astronomy Centre on low limit ( jParise et al. || 2003i HDO/H2O < 0.01 in ices), which 

behalf of the Particle Physics and Astronomy Research Council of the might be surprising if all species formed simultaneously on dust 

United Kingdom, the Netherlands Organization of Scientific Research, surfaces, 
and the National Research Council of Canada. 

** Member of the International Max Planck Research School (IMPRS) This low deuterium enrichment of water, if confirmed, 

for Astronomy and Astrophysics at the Universities of Bonn and would be a very valuable constraint for astrochemical models 

Cologne. that strive to explain the chemical processes involved in the 
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formation of water. Studying the abundance and distribution 
of HDO in low-mass prestellar envelopes is moreover an 
essential astrochemical goal, as it would provide the early 
conditions from where comets form. The HDO/H2O ratio has 
been measured in several comets and found to be ~3xlO 
dBalsiger et al.ll 19951: iBockelee-Morvan et al.|[T99l iMeier etall 



Il998h . Up to now, the HDO emission has been studied in detail 
only in one low-mass protostar (IRAS 16293-2422, Parise et al., 
l2005ah . A radiative transfer analysis of HDO emission has 
shown that the abundance profile of HDO can be constrained 
in low-mass protostellar envelopes by means of the analysis 
of several HDO lin es, spanning different energy conditions 
(iParise et al.l l2005bl) . Now that the Herschel telescope is 
delivering its first data, which allows us to get a better measure 
of the water abundance profile in protostars, the time is ripe 
to study water fractionation in these environments. Here we 
present the modeling of HDO observations toward a second 
solar-type Class protostar, NGC1333-IRAS2A (IRAS2A). 
In Section 2 we present the observations and first results, 
in Section 3 we present the radiative transfer modeling, and 
discuss the results in Section 4. We finally conclude in Section 5. 



2. Observations and results 

2.1. The source 

IRAS2A is a solar-type Class protostar located in the NGC 
13 33 molecular cloud . The adopted distance of IRAS2A is 235 
pc (iHirota et al.l2008l) . Deuterium fractionation of formaldehyde 
and m ethanol in the envelope has been studied by Parise et all 
(120061) . The deuterium enrichments (HDCO/H2CO ~ 0.17 and 
CH2DOH/CH3OH ~ 0.62) are higher than in IRAS 16293-2422 
jParise et al II2OOI . 

The H2O line emission in IRAS2A was recen tly observed 
with the Herschel telescope (iKristensen et alJl20Toh . 



2.2. Observations 

Observations were carried out at three different telescopes. 
Using the IRAM 30m telescope, we observed the 80, 225, 
and 241 GHz lines (see Figure [1) toward IRAS2A, at position 
Q'2000=03''28"'55'.60 and (52ooo=31''14"'37'.00, on 2004 Nov. 
26 under PWV ~ 2mm. The focus was checked on Uranus, 
and the local pointing on 0333-1-321, leading to an uncertainty 
in the pointing of less than 3" (rms). We used the wobbler 
switching mode, with a throw of 90". Typical system temper- 
atures varied in the range 115-150 K at 80 GHz, 240-3 10 K at 
225 GHz, and 375-525 K at 241 GHz. During a second IRAM 
run, we observed an offset position relative to the source, to 
check for a possible contribution of the outflow. We targeted the 
80GHz and 225GHz lines on 2005 Apr. 9, toward the position 
a2000 =03''29"'01''.00 and (52ooo=31"14"'20'. 

The 464 GHz observations were retrieved from the JCMT 
archive (project M04BN06). The lo,i-Oo,o line was observed to- 
ward the position of IRAS2A as well as toward the outflow on 
2004 Sep. 27, under an opacity T225ghz=0.05. 

The 893 GHz Une was observed using the CHAMP+ multi- 
pixel receiver on APEX in Sep. 2008. The central pixel was cen- 
tered on position q-2ooo=03''28"'55'.40 and 52000=31"14'"35'.00. 
Local pointing was checked in CO(6-5) on IK-Tau. We used the 
wobbler symmetric switching mode, with a throw of 120", re- 
sulting in OFF positions at 240" from the source. The total in- 
tegration time is 95 minutes. The CO(6-5) line, observed in par- 



allel on the 650 GHz receiver, is cl early detected with a pe ak 
of order of 3 K, and is presented in IKristensen et al.] (120101) . 

Table [T] lists the characteristics of the HDO observations. 
The unit of flux here was converted from to Tmb using the 
beam efficiencies indicated in the Table, which we took from the 
IRAIv£] and JCMlQ websites. The beam efliciency for CHAMP+ 
was measured in Oct. 2007 on Mars. 
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Fig. 1. HDO energy levels. The transitions that were observed in 
this work are marked with arrows (frequency in GHz). 



2.3. Results 

The spectra observed toward IRAS2A and the offset outflow po- 
sition are presented in Figure |2] The transitions from 80 GHz to 
464 GHz are detected on-source, while we only have an upper- 
limit for the second ground transition, at 893 GHz. No line is 
detected at the offset position (see Tableland Figure |2l). 



3. Modeling 

Because no emission is detected at the targeted outflow position, 
we use an envelope model to fit the observational results, exclud- 
ing any possible contribution of outflows or shocks. To analyze 
the HDO dat a we used the ID Monte C a rlo co de, RATRAN, 
developed by iHogerheiide & van der (120001) . The adopted 
density, opacity, and temperature profiles are derived with the 
ID DU STY code (Ivez ic & Elit zur 1997,) followi ng the proce- 
dure o f |j0rgensen elan 120021) dVildiz et alJl2010l) . In addition, 
we assumed the free-fall velocity field, v — (^'^'^' )'- and adopted 
the somewhat arbitrary values for the central mass, M, - 0.07 
and 0.15 Mq, consistent with the low mass inferred for the cen- 
tral object by iBrinch et al.l ( l2009l) . We will discuss this choice 
in Section 4.1. The emission was modeled in terms of a jump 
model, where the fractional abundances of deuterated water, rel- 
ative to H2, in the inner part of the source {T > 100 K, assumed 
evaporation temperature) and in the outer part {T < 100 K) are 
two free parameters, x™° and x™°, respectively. We performed 
anx^ analysis for x™° ranging from lxlO-^° to 2x10"^ and for 



' http://www.iram.es/IRAMES/telescope/telescopeSummary/ 
telescope_suminary.html 
^ http://www.jach.hawaii.edu/JCMT/spectraLline/Standards/beameff.html 
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Table 1. Parameters of the HDO observations 



Telescope 


Transition 


Frequency 


Eup 


Beam 


Seff 


rms* 


Resolution 


peak 


Av 








GHz 


K 


n 




mK 


km/s 


mK 


km/s 


Kkm/s 


NGC1333-IRAS2A 


IRAM 


1 1,0-1 1,1 


80.578 


46.8 


31.2 


0.78 


11.6 


0.15 


17.8 


3.9±1.1 


0.07±0.02 


IRAM 


3l,2-22,l 


225.897 


167.7 


11.1 


0.57 


24.7 


0.10 


113.5 


4.2±0.3 


0.50+0.03 


IRAM 


2i,i-2i,2 


241.561 


95.3 


10.4 


0.46 


30.8 


0.19 


98.6 


4.1±0.5 


0.43+0.05 


JCMT 


lo,i-Oo,o 


464.924 


22.3 


10.8 


0.44" 


125 


0.20 


305.8 


6.2±0.8 


2.00+0.23 


APEX 


li.i-Oo.o 


893.639 


42.9 


7.0 


0.35 


210 


0.25 








Outflow-position 


IRAM 


li,o-li,i 


80.578 


46.8 


31.2 


0.78 


11.5 


0.15 






<0.03" 


IRAM 


3l,2-22,l 


225.897 


167.7 


11.1 


0.57 


21.1 


0.10 






<0.04" 


JCMT 


lo,i-Oo,o 


464.924 


22.3 


10.8 


0.44" 


241 


0.20 






<0.72" 



"Values as observed on Jupiter during the second semester of 2004 with the RxW-C receiver. 
''T„it scale 

^3cr upper limit of integrated flux with assumption that Av is 5 km/s. 




10 20 
Vlsr (km/s) 

Fig. 2. Observed HDO spectra. The words in the lower left side 
of each spectrum indicate the observed position: source (on 
source) and outflow (offset-position). 



ranging f r om 9x 10 to 5x10 ^. Contrarily to the study 



by Parise et aP (l2005ah . who only modeled the integrated flux 
of HDO lines, we here intend to model as well the line veloc- 
ity profiles. In order to get a higher signal-to-noise ratio and a 
consistent weighting of each spectra, we smoothed the observed 
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Fig. 3. and contours (1,2 and 3cr) for the;^^^^^ for the 
different central star mass (M,). The black-solid and red-dashed 
fines indicate tfie contours of M, = 0.07 M© and M« - 0.15 
Mq. The symbols "-h" correspond to the best-fit model of each 
case. 



spectra to similar resolutions between 0.4 - 0.6 km/s. The;\f^ was 
then computed on the profile of these smoothed spectra. Here 

the definition of is x ''^"''''"'^'^"''^''""'''' , where the sum is over 
each channel in each spectrum. The parameters of the fits and 
their results are listed in Table [2] The cr within this analysis 
includes the statistical errors and uncertainties in flux calibration 
(~20%), but does not include any uncertainty in the adopt ed col- 
lisional rate coefficients used in the excitation calculation dGreenl 
Il989h . The black-solid and red-dashed lines in Fig. [3] present 
the contours dehmitating the Icr (68.3%), 2cr (95.4%) and 3cr 
(99.7%) confidence inte rvals. These contours are derived with 
the method described in iLampton et al.l (Il976h . who shows that 
the AS' = y~ - random variable follows a distribution with 
p variables, p being the number of parameters. Here p=2 (X,„ 
and Xo„,), so that these contours correspond to -x\^^^+'^-'i, 
X~ ■ H-6.17andy^. -nil. 8. Here we use Icr , 2 cr , and 3cr in 

mm ^ mm ' ' 

analogy with gaussian-distributed random variables. 

The obtained minimum is not far from 1, which indi- 
cates that the model assuming a simple collapsing envelope is a 
reasonable description of the data. The inner and outer fractional 
abundance is well constrained by our data at the 3cr level. The 
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Table 2. Fitting parameters and their results 



Central star mass" 


Doppler parameter b" 


Best fit HDO abundances 




HDO abundances confidence interval (3cr) 


(Mq) 


(km/s) 


inner region outer region 




inner region outer region 


0.07 
0.15 


2.0 
1.0 


8.0 xlO"** 7.0x10^'° 
8.6 X 10-** 5.5 X 10-'° 


1.1 
1.1 


6.6 X 10-'*-1.0 X 10"' 9 X 10-"-1.8 x lO^' 
6.8 X lO '^-l.l x 10 ' 9 X 10-"-1.5 x 10"' 



" Fixed parameters 

Table 3. Comparison of HDO fractional abundance between NGC1333-IRAS2A and IRAS 16293-2422 



Source 



best fit 



confidence interval (3cr) 



HDO 
H2O 



HDCOt 
H2CO 



CHoDOH f, 
CH3OH 



Inner outer (3cr) 



IRAS2A" 
IRAS 16293"' 



8 X 10^*' 
1 X 10-' 



7x 10" 
1.5 X 10 



6.6 X 10-**- 
7 X 10-**- 



1.0 X 10-' 
1.3 X 10-' 



9x 10-"~1.8x 10-' 
< l.Ox 10-' 



>o.or 

0.03 



07+°"'' 

' -0.06 

< 0.002 



17+° '2 
0.15±0.07 



62+°" 
37+°-^'* 

"•-"-0.19 



Orion KL 



0.02'' 



0.14' 



0.04s 



"The HDO fraction al abundance was derived with M=0.07 Mq. 

iParise et al.l l |200^ 

''The fractional abundances of water taken here are x|J-° <1 x 10-^ and x^^^ =1 x 10-^, as constrained bv lVissej et al. in prep, from the In - Ooo, 

202 - 111, and 2^ — 2q2 H'*0 lines and the upper-limit on the In - Ooo H'^^O line presented in Kristensen et al. (2010). 

' Parise et alj ( l2005a'l 

' Bergin et alTj201(]h 

^ Tumejjl99(T) 

' jjacQ etalJ ( li993h 




-10 10 20 



Vlsr (km/s) 

Fig. 4. Comparison of the observed spectra smoothed to 0.4 - 0.6 
km/s spectral resolution with the results of our best-fit models for 
two different masses of the central protostar The observed line 
emission is shown in black lines in each spectrum. The red and 
green lines show the results of the models with the assumptions 
that M, is 0.07 Mq and 0.15 Mq, respectively. 



smoothed spectra are overlaid with the modeling results of the 
best fit with two different central-star masses in Fig.|4] 



4. Discussion 

The first results of our model are that the observations are well 
described by the adopted envelope model, and that the HDO 
abundance profile displays a jump by two orders of magnitude 
between the cold and the warm envelope for the best fit. Below 
we first discuss the choice of the velocity field and the Doppler 
^-parameter, which is the turbulence broadening. 

4. 1 . The high turbulence broadening 

Although the minimum value shows that the simplest model 
(i.e. an infalling envelope) describes our data well, the required 
turbulence broadening b of 2 km/s may seem unrealistically high 
for the best model (M=0.07 M0). We note however that this 
b parameter is degenerate within the noise of our observations 
with the choice of the free-fall velocity profile in the envelope. 
Indeed, increasing the central mass (which is anyway mainly un- 
constrained) to 0.15 Mq allows us to reproduce the linewidths 
with a significantly lower b = I km/s. The minimum x^^^ in 
this case is the same as the model with central mass 0.07 Mq, 
and the x"°° and x™° of the best fit model are 8.6 x 10-*^ and 

m out 

5.5x10-'°, which are still within the Icr of previous results (Fig. 
|3]l. This shows that the constraints on and x™° do not de- 
pend much on the chosen velocity profile, and that the turbulence 
is barely constrained in the models. 

4.2. The fractional abundance of HDO 

The low value of x^^^ - 1-1 for the best-fit model shows that 
the assumption of a jump profile for the abundance in a spher- 
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ical envelope is appropriate to fit our data. The abundance of 
HDO jumps by at least a factor of 44 at the 3cr confidence 
level (see Fig. O. Ices are therefore shown to evaporate from 
the grains in the inner part of the envelope in IRAS2A. This 
prope rty has also been shown on HDO and H2O in IRAS 16293- 
2422 jParise et al.ll2005al: ICeccarelU et al.ll2000h . Moreover, the 
best-fit HDO abundances of the two cases are also similar to the 
values of HDO fractional abundances in IRAS 16293-2422 (see 
Table O. This suggests that the mechanisms to form HDO could 
be the same in both these solar-type class protostars, even if the 
two sources are located in different star-forming regions (Perseus 
v.s. Ophiuchus), and that the case of IRAS 16293-2422 is not 
peculiar. The present study provides the second detailed multi- 
transition HDO analysis toward a low-mass protostar, and the 
study of a greater number of sources would be essential to gen- 
eralize this result. 
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Fig. 5. Comparison of the HDO emission with H2O emission 
with/without wing. The red line indicates the emission line 
of the HDO transition at 225 GHz . The black and blue line 
show the Une emission of the H2O transition at 987 GHz with 
and without outflow component. The H2O data are taken from 
[Kristensen et al, (.201 ft) . 



4.3. Water deuterium fractionation in low-mass protostars 

In order to derive the HDO/H2O ratio in the IRAS2A envelope, 
we need to infer the H2O abundance profile in the quiescent en- 
velope. Figure |5] compares the spectra of the 225 HDO line and 
the 987 H2O Hne (Kristensen et al. 2010). Both the total wa- 
ter emission and the emission where the broad outflow Gaussian 
component was removed are shown. In this figure we multiplied 
the intensity of 225 HDO line with a factor of 2.2 to compare 
it with the H2O emission. Obviously the H2O line is broader 
than the HDO line (FWHM = 9.5 km s ' versus up to 6 km 
s'^ for HDO), even where the broad outflow component was 
removed. This medium-broad component has been interpreted 
to come from small-scale shocks because of the interaction of 
the jet and wind with the de nse inner envelope on scales of 
~ 1" jKristensen et alj r2010h . In addition, the H2O emission 
from the warm inner envelo pe was suggested to be optically 
thick (iKristensen et alJl2010l) . Therefore, the H2O lines do not 
probe the quiescent envelope and HDO lines may not probe the 
same gas. 



To obtain upper limits on the quiescent H2O abundance, 
Visser et al. in prep, used deep integrations on various H'**0 lines 
presented by Kristensen et al. (2010). The non-detection of nar- 
row H'**0 features limits the inner H2O abundance to < 10"^, 
whereas the presence of narrow deep H2O self-absorption gives 
an outer H2O abundance of ~ 10"^. 

The values of the HDO/H2O ratio within the present con- 
straints are listed in Table [3] The HDO/H2O ratios in the inner 
warm and the outer cold envelope are found to be more than 1 % 
and are in the range of 0.9%-18% (3cr). The D/H ratio of water 
in the outer envelope is roughly one order of magnitude lower 
than the ratio of methanol, which confirms the trend measured 
in IRAS 16293-2422. Obviously, better constraints on the H2O 
abundances in the inner warm envelope are required to see if 
the deuterium enrichment in water is significantly lower than in 
methanol. On the other hand, the similar HDO/H2O ratio and the 
one order magnitude higher CH2DOH/CH3OH ratio of IRAS2A 
and Orion KL suggest that the environment plays a role in the 
chemistry. 

The HDO/H2O ratio in the inner envelope is found to be 
more than one order of magnitude higher than the ratio measured 
in comets, as for IRAS 16293-2422. If this high deuterium en- 
richment of water is typical of low-mass protostellar envelopes, 
and if the HDO/H2O ratio measured in the cometary coma is 
representative of the cometary ice composition, then this would 
imply that the deuterium fractionation of water is reprocessed at 
some point between the protostellar envelope and the cometary 
ice. Reprocessing of the gas in the hot corino warm environ- 
ment is expected to reduce the deuterium fractionation on typi- 
cal timescales of 10^ yrs in environments of density ~10'' cm"^ 
(ICharnlev et al.l[T997l) . However, in this case, re-condensation 
of the gas on the grains would be required. Isotopologue ex- 
changes are also pos sible at the surface of the grains (see e.g. 
iRataiczak et"ani2009l) . without requiring evaporation. This may 
however not be enough to deplete the deuterium level in water, 
because water is the main component of water ices. A mecha- 
nism to get rid of deuterons from the grain surface would need to 
be invoked. Because of the high abundance of water, this mech- 
anism should involve efficient exchange reactions. A suggestion 
would be to study in the laboratory the possibility of deuteron 
exchange between HDO and H2 at the grain surfaces. 

Improving the observational constraints on the HDO/H2O ra- 
tio in protostellar envelopes is essential to allow us to understand 
the implications in terms of chemical evolution in the protosolar 
nebula. Figure |6] shows the RATRAN modeling prediction for 
two H^^O lines, for a set of reasonable water abundances. We 
made the assumption that the standard ortho/para ratio is 3 and 
the i**0/'*'0 ratio is 2.05xl0"\ Figure |6h shows that the H^'^O 
1 189 GHz line emission is much more sensitive to the inner frac- 
tional abundance of water than the Hj'^O 1101 GHz line. This 
line is therefore a good candidate to target with the Herschel 
telescope to further constrain the inner water abundance. 

5. Conclusion 

We presented the observation of five HDO lines toward the solar- 
type class protostar NGC1333-IRAS2A using the IRAM 30m, 
JCMT and APEX telescope. Four of them are clearly detected, 
while the second ground-state transition line at 893 GHz is not. 
Because there are no detections at the outflow position at 80, 
225, and 464 GHz, we assumed for the modeling that most emis- 
sion comes from the envelope. 

We modeled the emission with the RATRAN radiative trans- 
fer code, assuming a step profile for the HDO abundance. We 
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Fig. 6. Model predictions for the chosen ^[2*^0 lines. X,„ and Xo„, 
represent the water (H2O) abundance in the inner warm (> 100 
K) and outer cold ( < 100 K) envelope. A standard ortho/para 
ratio of three was assumed. 

derived the HDO abundance in the inner and outer parts of the 
envelope to be x™° = 8 x lO"*^ and x™° = 7 x lO"'". This 
result shows that the HDO abundance is enhanced in the inner 
envelope because of the ices' evaporation from the grains, as for 
IRAS 16293-2422. Moreover, the values are similar to those in 
IRAS 16293-2422, which suggests that the pathway of the for- 
mation of HDO is the same in low-mass class protostars. The 
study of a larger sample of sources is needed to generalize this 
result. 

The HDO/H2O abundance ratio is found to be > 1 % in the 
inner envelope. A better constraint on the H2O abundance in the 
inner envelope is required to derive the deuterium enrichment in 
the warm envelope. 
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